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ABSTRACT 

Dynamical considerations, presented herein via analytic scalings and numeri- 
cal experiments, imply that Earth-mass planets accreting in regions that become 
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habitable zones of M dwarf stars form within several million years. Temperatures 
^ ■ in these regions during planetary accretion are higher than those encountered by 

the material that formed the Earth. Collision velocities during and after the 
prime accretionary epoch are larger than for Earth. These factors suggest that 
^ ■ planets orbiting low mass main sequence stars are likely to be either too dis- 
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tant (and thus too cold) for carbon/water based life on their surfaces or have 
abundances of the required volatiles that are substantially less than on Earth. 

Subject headings: (stars:) planetary systems: formation; planets and satellites: 
formation; astrobiology 



1. Introduction 

M (and early L) dwarfs are the lowest mass stars, and the smallest and least luminous 
members of the stellar main sequence. They are also by far the most numerous class of 
stars in our galaxy, with ~ 75% of the stars in the extended solar neighborhood being M 



dwarfs (IReid et al.ll2002l ). Thus, if planets orbiting M stars can be habitable, the number 
of potentially inhabited worlds within our galaxy could be much larger than if only planets 
orbiting solar type stars are capable of hosting life. The presence of liquid water on a planet 
for long periods of time is considered to be a requirement for the origin and evolution of life 
as we know it. It is shown herein that planets within M star habitable zones (HZs) formed 
rapidly while their star was quite luminous and are probably water-deficient. Thus, the 
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dynamics of planetary accretion, together with the physics of stellar evolution, cast doubt 
on the suitability of planets orbiting M stars to host life. 

Stars on the main sequence have relatively long periods of very slowly varying luminosity, 
potentially providing stable radiation fluxes for the origin and evolution of life. Giant planets, 
brown dwarfs, and post-main sequence stars vary substantially in luminosity, and thus are 
far le ss conducive energy sources for life on the surface of a planetary body (cf. lAdams et al. 
20051 . for a possible exception involving a class of low-mass post-main sequence stars in the 
very distant future). 

The main sequence luminosity of stars varies roughly as Mf , where M* is stellar mass. 
Thus, the locations of HZs (defined to be orbits in which liquid water can be present on the 
surface of an Earth- like planet) are much cl oser to low mass main sequence stars than they 
are to high mass stars feasting et al.lll993l ). As stellar nuclear fuel varies in proportion to 
M*, stellar lifetime goes as M~ 3 , so stars more than about twice as massive as the Sun are 
probably not sufficiently long-lived to provide the energy for the development of advanced 
life. 

The lowest m ass stars are very long-lived, but present a different set of challenges to life 
(jScalo et al.ll2007l ). The most studied of these problems is that the proximity of HZs to very 
low mass stars implies that the rotation of a HZ Earth-like planet would be tidally altered 
to the sy nchronous state, wit h one hemisphere permanently illuminated and the other in 
the dark (jKasting et al.lll993l ). A thin atmosphere of such a planet would freeze out onto 
the planet's night hemisphere, but a sufficiently massive atmosphere could transport enough 
heat to the dark side to p revent atmospher ic collapse, thus allowing habitable regions on 
the planet's lit hemisphere (IJoshi et al .11199 71 ) . Perturbations from another planet could also 
maintain a HZ planet's orbital eccentricity against tidal locking, especially if the HZ planet 
lacked a substantial permanent deformation (and thus could rotate just slightly faster than 
synchronous, in equilibrium with the tidal forcing experienced by an eccentric planet). 

However, M dwarf stars present other challenges to planetary habitability. As shown 
below, a particularly severe problem is the difficulty that a 1 M e (Earth mass) planet in the 
HZ of an M dwarf has in accumulating and retaining water and other volatiles. The Earth 
itself is volatile-poor. Oceans and other near-surface reservoirs of H2O comprise less than 
0.03% of our planet's mass, with the mantle containing a comparable amount of water (K. 
Zahnle, personal communication, 2006). Thus, Earth's total water abundance is very small 
compared to the cosmic rock:H 2 ratio of roughly unity. Earth accreted from planetesimals 
that condensed over a wide range of heliocentric distances (jWetheriillll994l ) . The temperature 
distributions of protoplanetary disks are not tightly constrained. Thus, the best estimates 
of the volatile contents of planetesimals as a function of heliocentric condensation distance 
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are provided by analysis of meteorites that are analogous to asteroids of known heliocen- 
tric distance and by m atching; Earth's volatile inventory to planetary accumulation models 
( IMorbidelli et al.ll2000l ). Primitive meteorites that emanate from the inner asteroid belt (just 
beyond 2 AU) have water abundances < 0.1%, whereas those from more distant regions of 
the asteroid belt can have more than 100 times this much water. Accretion simulations 
imply that a nontrivial amount of material from the region of the main asteroid belt reached 
Earth, and that m ost of Earth's water probably came from p lanetesimals which condensed 
beyond ~ 2.5 AU (iLunine et al.ll2003l ; iRaymond et al.l I2004J ) . Although current planetary 
accumulation models cannot be viewed as providing ab initio estimates of the accretion of 
water by terrestrial planets, models of this type are well-suited for comparing the relative 
amounts of volatiles accreted by terrestrial planets around stars of differing masses. 

Stars are substantially more luminous during their formation epoch than they are once 
they reach the main sequence. Low mass stars are significantly more luminous than t heir 
main sequence luminosities for longer than are solar mass stars (IStahler fc Pallal 120051 ). so 
ultimately habitable regions around M stars are hotter at young stellar ages than are regions 
that are to be subjected to the same stellar flux around sunlike stars. Additionally, energy 
diffusion is slower in inner protoplanetary disks because of the greater optical thickness, so 
there is likely to be an even larger excess of temperature within young circumstellar disks 
near smaller mass protostars. The loca tion of the snow line within protoplanetary disks 



orbiting small stars ([Kennedy et al.ll2006l ) is thus more distant in proportion to the eventual 
location of the HZ than is the case for solar mass stars. 

The timescales of accretion of planets of a given mass in the HZs of M stars are shorter 
than those of comparable planets around more massive stars. This is because orbital periods 
are shorter, planetesimals are closer to one another and they occupy a larger fraction of their 
Hill spheres. Moreover, orbital velocities are faster, implying higher impact speeds, so late 
accretion of volatile-rich bodies that condensed farther from the star may well remove more 
atmospheric gasses and water than they provide. Thus, while lower mass pre-main sequence 
stars remain more luminous than their main sequence luminosities for a longer time, Earth- 
mass planets forming in HZs around these stars accrete more rapidly, of material that is 
likely less water-rich, and within a dynamical environment in which they are more likely to 
lose atmospheric volatiles via impact erosion. The dynamical characteristics of terrestrial 
planet accumulation around M stars are quantified in Section 2. 
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2. Accretion Timescales & Impact Velocities 



For stars of solar mass and smaller, main sequence luminosity, L*, varies with stellar 
mass, M*, roughly as: 



M 



where the symbol ^ signifies that the proportionality relationship is only approximate^. 
Planetary temperature is proportional to (L*/r 2 ) 1//4 , where is the distance from the star. 



The distance of the HZ from the star, Tjjz, thus varies as: 

r HZ KLl' 2 ^Ml 



(2) 



Terrestrial p lanet accretion proceeds until a stable configuration is reached (jLissauerlll995 



Laskarl 2000), and the separation of planets required for stability varies as [M v j 'M*) 2 / 7 
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19801 ) . Thus, the surface mass density of the ensemble of solids within the proto- 
planetary disk, er, required for local accretion of planets of a given mass (e.g., that of Earth) 
in the zone that becomes habitable varies as 
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When expressed in terms of surface density in the HZ, the value of a necessary for the 
formation of earthlike planets is, therefore, a very ste ep function of stellar m ass. But standard 



minimum mass models of our protoplanetary dislcl ( Weidenschillingl 1 1 9 771 ) give a 
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which is only slightly less steep than the r^^ 7 dependence in Eq. (3). Note that masses 



of protoplanetary disks around low mass stars are poorly constrained , but observations 
suggest that they ty pically are lower than those around 1 M stars (IScholz et al.l 12006 



Muzerolle et al.l 120061 ) . If the surface densities of disks in ultimately HZs around small stars 
do not have substantially higher surface mass density than the minimum mass proto-solar 



1 The actual slope of the mass-luminosity relationship for 0.5 - 1 stars is slightly ste eper, between 
and M+, but it flattens out to vary ~ Mj for lower masses ([Hillenbrand fc White! 12004 , and J. Scalo, 



personal communication, 2006). Therefore, the HZ is a little closer to K and early M stars (M* > 4 M Q ), 



than predicted using Eq. (1), but a somewhat farther from late M stars (M* < 3 
scaling. 



1 M ) 



than given by this 



2 These models are not constrained by data from regions as close to the Sun as the HZ distance is from 
an M star, because our Sun has no planets that orbit so near it. Also, flatter surface density models are 



derived in other minimu m mass solar n ebula models (jDayisJ [2005J) , and such flat profiles better fit giant 
planet formation models ( Lissauerlll987| ) as well as most evolutionary models of protoplanetary disks. 



- 5 - 



disk at 1 AU, then planets in the HZs of M stars are small unless they accreted large amounts 
of material that migrated inwards from the outer disk or they themselves migrated inwards 
substantially. Low-mass planets have shallower gravitational potential wells, and thus are 
more susceptible to atmospheric loss than are more massive planets. 

The growth rate of a terrestrial planet is proportional to the surface density of the disk 
and the orbital frequency; ot her factors rela ted to gravitational enhancement in accretion 
cross-section also play a role (jSafronovl Il969l ) . but this enhancement is of order unity dur- 
i ng the final high-velocity phase o f planetary growth, which dominates accretion timescales 
( lLissauerlll987l ; lAgnor et al.lll999l ). Growth times thus scale as 
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According to Eq. (4), Earth-mass planets that accrete within HZs of 0.5 M stars should 
form in several million years; around 0.25 M Q stars, the process should require < 10 5 years. 
Migration inwards requires the presence of a massive disk. Optically thick disks around 
young stars have a broad range of lifetimes, alt hough most d i sappe a r 1-5 Myr after t he 
stellar photosphere becomes visible (Figure 5 of iBriceno et al.l (120071 ); iMeyer et al.l (120071 )). 
Thus, whether it has grown in situ or migrated inwards from a greater distance, the planet 
should be in or very near its final orbit within the HZ within 10 7 years of the star's formation. 

A slightly less steep scaling with stellar mass than given by Eq. (4) is probably more 
appropriate to account for two physical effects: (1) The amount of gravitational focus- 
i ng is reduced for accreting bodies which occupy a greater fraction of their Hill spheres 
( jGreenzweig fc Lissauerl 1 19901 ) . (2) The formulae used to derive Eq. (4) assume crossing 
orbits, whereas the very late stages of terrestrial planet growth are dominated by the time 
required for chaotic perturbations to excite sufficient eccentricities for continuation of orbital 
crossing. 

The growth of terrestrial planets within the HZ of a | star (which has a mai n 



sequence luminosity very close to (l/3)~ 4 that of a 1 M Q star ([Hillenbrand fc White! 120041 )) 
was modeled numerically to test the scaling derived above. The initial disk was analogous 
to that used for simulatio ns of terrestrial planet grow th around a 1 M Q star as well as 



within binary star systems (jOuintana fc Lissauerl 12006k and references therein). As in these 
previous studies, the simulations began with 140 'planetesimals' of mass 0.00933 M ffi and 
14 'planetary embryos' of mass 0.0933 M e , and the density of each body was taken to be 
3 g/cm 3 . In order to leave the illumination of bodies that were at 1 AU in the initial disk 
unchanged according to Eq. (2) and maintain expected planetary sizes as per Eq. (3), the 
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initial semimajor axes of these 154 bodies were scaled as 

J3 2/7 ) 



9 



(5) 



where a p are the values used in the previous simulations referenced above and both a and a p 
are in AU. The other 5 initial orbital elements of each of the 154 accreting bodies were as in 
previous simulations. Bodies were removed from the integrations if they approached closer 
to the star than 0.01 AU, or if their distance from the star exceeded 10 AU. The integrations 
were performed using the hybrid symplectic integrator within the Mercury integration pack- 



age (IChambersI Il999l ) . The timestep was taken to be one hour. Deterministic chaos implies 
that integrations of this sort are valid in only in a statistical sense, so multiple numerical 
experiments were run. A pair of simulations, differing only by moving the initial position 
of one of the smaller bodies along its orbit by one meter, was performed without any giant 
planets, and four analogous runs were performed with giant planets similar to Jupiter and 
Saturn but with masses reduced by a factor of three and semimajor axes reduced by a factor 
of nine. 

Results of the four numerical simulations of planetary growth in the HZ of a | M star 
that has a pair of giant planets are displayed in Figure 1. The final systems have similar 
numbers of terrestrial planets as those formed in simulations of planetary growth around 1 
M Q stars on which the initial disk conditions were based. The final collision between two 
embryo-dominated bodies or loss of one such body occurred at 0.39 Myr, 0.38 Myr, 1.21 Myr 
and 0.67 Myr in these 2 Myr simulations; in comparison, the last reduction in number of 
embryos in t he 31 simulations of te r restri al planet growth (each for a time span of 200 Myr) 



tabulated in iQuintana &: Lissauerl (120061 ) occurred at times ranging from 18.7 Myr to 134 
Myr, with a median final loss time of 78 Myr. Growth rates around the | M star thus are 
more than 100 times as rapid as around the 1 M star, which implies a somewhat less steep 
scaling than Eq. (4), consistent with the above discussion. The simulations lacking giant 
planets resulted in more final terrestrial planets extending to greater astero centric distance, 



as in the case of a single star lacking in giant planet or stellar companions (IQuintana et al. 



20021 ). In these runs, the final collision between two embryo-dominated bodies or loss of one 
such body occurred at 0.11 Myr and 0.18 Myr. Note that ~ M Q stars are more than ten 
times as luminous during thei r first ~ 4 Myr as they are when they reach the main sequence 



(ID'Antona fc Mazzitellilll997h 



Planets grow by accretionary impacts. Such impacts can remove previously accreted 
volatiles that reside in a planet's atmosphere. Impact erosion is a favo red mechanism for 



expla ining the thinness of Mars's atmosphere relative to that of Earth (IMelosh fc Vickery 



19891 ). Orbital speeds, and thus impact velocities, are higher in HZs around low- mass stars 
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than they are at Earth, and such high collisional velocities may lead to more erosive impacts. 
Earth is believed to have accumulated most of its water and other volatiles from solid bodies 
that condensed outwards of the orbit of Mars and were ac creted by Earth toward s the end of 
or subsequent to the major phases of our planet's growth (IMorbidelli et al.ll2000l ). Accretion 
from a similarly-scaled distance around other stars would result in approach velocities, v^, 
that vary in proportion with the orbital velocity!^, 
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ocMT 1/2 . 



(6) 



The larger ratio of the distance to the ice line to rnz for M stars implies that fewer ice-rich 
planetesimals reach the HZ, and those that do have an even higher characteristic velocity 
than given by Eq. (6 ). The specific en e rgy o f an impact varies as the square of impact 
velocity. According to iMelosh fc Vickeryl (119891 ). substantial impact erosion of atmospheres 
requires collision velocities faster than twice the planet's escape speed. Parabolic comets 
with prograde orbits in the ecliptic plane impact Earth at < 2v e . However, all planetesimals 
with large semimajor axes impact a 1 M e planet orbiting at | AU from a | M Q star at > 2v e . 
Therefore, impacts in the HZs around small stars can significantly erode the atmospheres of 
< 1 M e planets. 



3. Discussion 

In sum, under nominal circumstances, planets in main sequence habitable zones around 
M stars are likely to be fully formed and in their final orbits by the time the gaseous circum- 
stellar disk has dissipated or several million years after planetesimal formation, whichever 
is later. If growth is in situ, dynamical and thermal factors imply that the planets are un- 
likely to have large volatile inventories, and planetary masses are likely to be small. The 
large collision speeds of impacting comets, as well as the high activity and luminosities of 
young M stars, may lead to substantial mass loss from planetary atmospheres, depleting any 
reservoirs of volatiles that planets within the HZs are able to accrete. 

So are M dwarf stars totally unsuitable hosts for the development of advanced life? 



3 The situation compares even less favorably for 'local' impacts during the late stages of the primary 
growth of planets that reach the same masses, as they must be separated by a larger fraction of their orbital 

2/7 

distance for stability about the smaller star. Thus, an additional factor of M* must be included to 
account for the greater eccentricities of planetesimals which are cleared from the entire region between two 
growing planets, yielding Voo JJM* n / 14 . 
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No, there are various ways out of the difficulties pointed out above. Theoretical estimates 
of planetary volatile inventories obtained from planet formation studies are not from first 
principles, but rather rely on normalization to the planets within o ur Solar System; a s 
accretion models suggest large stochastic variations in volatile delivery (jO'Brien et al.ll2006l ). 
it is possible that the terrestrial planets within our Solar System are near the volatile-poor 
end of this distribution, and thus th e normalization used herein is incorrect. A water-rich 
planet formed farther from the star ([Kennedy et al.ll2006l ) could migrate inwards to the HZ 
while the gaseous disk is still present; such a planet could initially have enough water that 
oceans could be retained even if it suffered significant losses during the star's young active 
phase. A planet could be placed on an eccentric orbit that circularizes on the timescale 
it takes the star to reach the main sequence or somewhat longer, or be scattered inwards 
by a fellow planet after the star reaches main sequence and then be circularized by stellar 
tides. But such scenarios require precisely the right amount of initial water or just the right 
dynamics. Thus, while it is likely that some of the hundreds of billions of M dwarf stars in 
our galaxy have planets with temperatures, masses and compositions similar to Earth, the 
number of such planets is probably small, and Sun-like stars, despite being considerably less 
numerous, may well be the hosts of far more habitable planets. 
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Fig. 1. — The temporal evolution of four protoplanetary disks incorporating the ultimately 
habitable zone around a | M Q star. Two giant planets, one-third the masses of Jupiter 
and Saturn, respectively, are also included. The planetary embryos and planetesimals are 
represented by circles whose sizes are proportional to the physical sizes of the bodies. The 
evolution of one disk is shown in black at 12 different times. The green and red filled circles 
and blue open circles show the results at late times of the three runs with one planetesimal 
initially moved by 1 - 3 m along its orbit. The horizontal locations of the circles show 
the orbital semimajor axes of the bodies in the disk in AU, and vertical positions plot 
their eccentricities. The initially dynamically cold disk heats up during the first 5 x 10 4 
yr, especially in the outer region, where the perturbations of the giant planets are the 
greatest. By 4 x 10 5 yr into the simulation, 3-5 terrestrial planets have formed, with 
0-2 planetesimals remaining. After 2 x 10 6 yr, each system contains 2-4 terrestrial 
planets. Note that the HZ is around 0.11 AU and that the largest planets formed have 
masses comparable to that of the Earth. 



